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This  article  aims  to  develop  a  thermodynamic  modelling  and  optimization  framework  for  a  thorough 
understanding  of  the  optimal  integration  of  fuel  cell,  gas  turbine  and  other  components  in  an  ambient 
pressure  SOFC-GT  hybrid  power  plant.  This  method  is  based  on  the  coupling  of  a  syngas-fed  SOFC  model 
and  an  associated  irreversible  GT  model,  with  an  optimization  algorithm  developed  using  MATLAB  to 
efficiently  explore  the  range  of  possible  operating  conditions.  Energy  and  entropy  balance  analysis  has 
been  carried  out  for  the  entire  system  to  observe  the  irreversibility  distribution  within  the  plant  and  the 
contribution  of  different  components.  Based  on  the  methodology  developed,  a  comprehensive  parametric 
analysis  has  been  performed  to  explore  the  optimum  system  behavior,  and  predict  the  sensitivity  of 
system  performance  to  the  variations  in  major  design  and  operating  parameters.  The  current  density, 
operating  temperature,  fuel  utilization  and  temperature  gradient  of  the  fuel  cell,  as  well  as  the  isentropic 
efficiencies  and  temperature  ratio  of  the  gas  turbine  cycle,  together  with  three  parameters  related  to 
the  heat  transfer  between  subsystems  are  all  set  to  be  controllable  variables.  Other  factors  affecting  the 
hybrid  efficiency  have  been  further  simulated  and  analysed.  The  model  developed  is  able  to  predict  the 
performance  characteristics  of  a  wide  range  of  hybrid  systems  potentially  sizing  from  2000  to  2500  W  m-2 
with  efficiencies  varying  between  50%  and  60%.  The  analysis  enables  us  to  identify  the  system  design 
tradeoffs,  and  therefore  to  determine  better  integration  strategies  for  advanced  SOFC-GT  systems. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

To  overcome  the  threats  posed  by  climate  change  and  energy 
security,  advanced  clean  energy  technologies  are  urgently  needed. 
Hybrid  fuel  cell  systems  have  emerged  as  attractive  power  genera¬ 
tion  solutions  with  great  promise  for  high  energy/power  efficiency 
with  reduced  environmental  emissions,  allowing  the  shift  from 
a  fossil  fuel-based  economy  to  a  new  paradigm  in  a  progressive 
manner  [1,2].  Amongst  the  various  types  of  fuel  cell  hybrids,  the 
combination  of  solid  oxide  fuel  cell  (SOFC)  and  gas  turbine  (GT)  has 
been  identified  as  a  key  technology  superior  to  many  other  options 
due  to  its  fuel  flexibility  and  ultra-high  energy  conversion  efficiency 
[1-3].  To  date  several  hybrid  configurations  have  been  presented 
or  patented  [4-6],  along  with  theoretical  studies  which  involve 
the  thermodynamic  analysis,  design  and  performance  modelling 
of  SOFC-GT  hybrids  through  individual  and  collaborative  efforts 
[7-17]. 

In  spite  of  the  advantages  of  the  SOFC-GT  hybrid  system,  many 
technical  barriers  have  to  be  overcome  for  the  successful  com¬ 
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mercial  development  of  this  power  generation  technology.  These 
barriers  are  particularly  evident  in  the  smooth  integration  of  the 
fuel  cell  and  GT,  since  the  GT  has  a  strong  coupling  with  the  top¬ 
ping  SOFC  [18].  Therefore,  design  parameters  need  to  be  varied 
to  determine  performance  sensitivity,  and  fundamental  analysis  is 
still  needed  to  determine  the  practical  operating  constraints  and 
the  conditions  for  most  efficient  operation.  Although  there  have 
already  been  a  number  of  papers  published  in  the  area  of  hybrid 
SOFC-GT  system  analysis,  most  of  them  are  generally  focused  on 
part-load  or  design  point  performance  simulations  instead  of  opti¬ 
mization.  No  optimization  strategy  has  been  proposed  so  far  from  a 
systematic  integration  perspective.  Moreover,  all  the  research  for¬ 
merly  done  was  limited  to  either  a  particular  system  layout,  or  a 
particular  operating  or  off-design  condition. 

System  modelling  combined  with  thermodynamic  optimiza¬ 
tion  can  be  a  valuable  tool  in  technological  research,  providing 
indications  of  technical  feasibility,  identifying  ways  to  improve 
efficiency,  and  determining  the  best  configuration  and  conditions 
for  an  integrated  power  plant  [19].  To  make  a  thorough  investi¬ 
gation  of  the  optimal  integration  amongst  a  fuel  cell,  gas  turbine 
and  other  system  components,  a  generic  thermodynamic  mod¬ 
elling  and  optimization  framework  is  developed  in  the  present 
study  to  systematically  optimize  the  ambient  pressure  SOFC-GT 
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hybrid  cycle  as  an  important  follow-up  research  effort  based  upon 
previous  system  modelling  work.  This  methodology  is  based  on 
the  coupling  of  a  coal  syngas-fed  SOFC  model  and  an  irreversible 
gas  turbine  cycle  model,  with  an  optimization  algorithm  devel¬ 
oped  using  MATLAB  to  efficiently  explore  the  range  of  possible 
operating  conditions.  In  order  to  observe  the  irreversibility  dis¬ 
tribution  within  the  plant  and  the  contribution  of  different  cycle 
components,  a  comprehensive  energy  and  entropy  balance  analy¬ 
sis  has  been  conducted  for  the  entire  system.  Multi-irreversibilities 
coming  from  the  electrochemical  polarization,  parasitic  power  con¬ 
sumption,  finite-rate  heat  transfer,  irreversible  compression  and 
expansion,  and  heat  loss  to  the  environment  are  all  taken  into 
account.  The  optimal  efficiency  and  power  density  of  the  hybrid 
system  have  been  deduced  as  a  function  of  various  system  param¬ 
eters.  The  current  density,  operating  temperature,  fuel  utilization 
factor  and  temperature  gradient  of  the  fuel  cell,  as  well  as  the  isen- 
tropic  efficiencies  and  temperature  ratio  of  the  gas  turbine,  together 
with  two  parameters  related  to  the  heat  transfer  between  subsys¬ 
tems  are  all  set  to  be  controllable  variables.  Detailed  numerical 
simulations  have  been  performed  to  explore  the  optimum  sys¬ 
tem  behavior  and  parametric  sensitivity  analysis.  It  is  important  to 
notice  that  the  approach  presented  here  is  based  on  an  integrated 
multi-scale  thermodynamic  modelling  and  optimal  design  concept, 
which  provides  a  new  perspective  to  understand  the  hybrid  fuel  cell 
system  with  reduced  complexity  and  enhanced  generality.  There¬ 
fore  one  of  the  most  useful  features  of  the  methodologies  is  the 
ability  to  use  the  model  to  systematically  evaluate  the  performance 
potential  for  a  wide  variety  of  configurations  as  well  as  the  sys¬ 
tem  design  tradeoffs.  This,  in  turn,  assists  the  optimization  and 
development  of  new  energy  technologies. 

2.  System  configuration  and  description 

The  hybrid  coupling  derived  from  the  integration  of  high 
temperature  SOFC  and  gas  turbine  technology  has  been  studied 
extensively  [20],  however,  the  prediction  of  the  overall  response 
of  the  system  performance  to  the  larger  number  of  operating  con¬ 
ditions  remains  a  rather  complex  task.  To  develop  an  effective 
optimization  strategy,  it  is  important  to  understand  the  system 
design,  component  interactions  and  operating  requirements. 

2.1.  System  configuration 

A  fuel  cell  and  a  gas  turbine  may  be  configured  in  several  differ¬ 
ent  ways.  They  can  be  built  to  operate  via  either  direct  or  indirect 
thermal  coupling.  Studies  indicate  that  atmospheric  system  config¬ 
urations,  in  which  the  thermal  energy  is  utilized  indirectly  through 
a  heat  exchanger  allowing  the  selection  of  GT  pressure  to  be  inde¬ 
pendent  of  the  fuel  cell  pressure,  would  be  less  complex  to  develop 
and  quicker  to  implement  than  pressurized  ones.  Since  the  fuel  cell 
and  gas  turbine  operate  with  different  working  fluids  and  at  dif¬ 
ferent  pressures,  the  ambient  pressure  systems  have  the  potential 
to  be  less  expensive  and  could  accommodate  a  wider  variety  of  gas 
turbines  as  well  as  fuel  cells  [21  ].  For  these  reasons,  an  atmospheric 
system  layout  has  been  chosen  to  be  examined  in  the  present  study. 
The  basic  concept  of  this  SOFC-GT  hybrid  is  depicted  in  Fig.  1, 
where  the  fuel  cell  acts  as  the  core  system  component  and  gen¬ 
erates  most  of  the  system  power.  Air  is  supplied  by  a  blower  and 
preheated  prior  to  entry  to  the  SOFC.  The  product  gas  of  SOFC  is 
sent  into  afterburner,  where  the  un-reacted  fuel  is  burnt  with  part 
of  the  excess  air.  Downstream  of  the  fuel  cell,  the  gas  turbine  is 
modelled  as  a  compressor  and  a  turbine  mechanically  linked  via 
a  common  shaft,  and  thermodynamically  coupled  to  the  fuel  cell 
via  a  primary  heat  exchanger  (HEX1 ).  The  air  leaving  the  compres¬ 
sor  is  heated  up  by  the  high-temperature  exhaust  gas  of  the  fuel 


Fig.  1.  The  schematic  diagram  of  an  atmospheric  SOFC-GT  hybrid  system. 


cell  which  carries  a  significant  amount  of  thermal  energy.  This  hot 
air  enters  directly  into  the  turbine  and  expands  to  produce  power. 
As  the  most  noticeable  difference  in  configurations,  two  supple¬ 
mentary  heat  exchangers  (FIEX2  and  FIEX3)  are  adopted  to  serve 
as  pre-heaters,  recovering  heat  from  both  the  exhaust  gas  of  FIEX1 
and  the  turbine  outlet  for  preheating  the  inlet  air  and  fuel  before 
they  enter  the  fuel  cell.  The  final  exhaust  gas  rich  in  carbon  diox¬ 
ide  from  the  hybrid  system  can  be  made  capture  ready  after  drying 
(water  removal)  and  compression. 


2.2.  Model  assumptions 

Mathematical  modelling  is  an  essential  tool  for  the  analysis  and 
design  of  energy  conversion  systems.  This  is  particularly  appro¬ 
priate  for  the  SOFC-GT  hybrids,  where  localized  experimental 
measurements  are  difficult  due  to  the  high  operating  temperature. 
To  develop  a  mathematical  model  representing  the  electrochemical 
and  thermodynamic  characteristics  of  the  entire  system,  a  number 
of  simplifications  and  assumptions  are  made  to  enable  the  analysis 
[7,19,21-30]: 


•  Both  the  fuel  cell  and  gas  turbine  are  assumed  to  be  operated 
under  steady-state  conditions. 

•  Fuel  cell  reactions  are  assumed  to  be  in  equilibrium. 

•  Syngas  consists  of  the  following  gas  species,  j  =  {H2,  CO,  CH4,  C02, 
H20,N2}. 

•  Air  that  enters  the  fuel  cell  consists  of  79%  N2  and  21%  02. 

•  The  cathode  and  anode  inlet  temperature  of  the  fuel  cell  are 
assumed  to  be  equal. 

•  The  cathode  and  anode  exit  temperature  of  the  fuel  cell  are 
assumed  to  be  equal. 

•  There  is  a  temperature  gradient  (AT)  across  the  SOFC.  The  tem¬ 
perature  of  the  solid  structure  (T)  is  homogeneous  and  midway 
between  the  inlet  and  exit  temperatures. 

•  All  gases  behave  as  ideal  gases. 

•  Gas  leakage  is  negligible. 

•  Heat  loss  to  the  environment  occurs  only  in  the  fuel  cell. 


With  the  help  of  these  assumptions,  the  system  model  will  be 
constructed  and  the  governing  equations  representing  all  modelled 
components  of  the  hybrid  cycle  will  be  given  in  the  following  sec¬ 
tions.  Each  of  the  system  components  are  modelled  individually 
and  integrated  to  form  the  overall  hybrid  system. 
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3.  Model  formulation 

3  A.  Solid  oxide  fuel  cell  fuelled  with  coal  syngas 

SOFCs  are  devices  that  continuously  convert  chemical  energy 
directly  into  electrical  energy  through  electrochemical  reactions. 
Various  fuel  options  are  feasible  for  SOFC  operation  [1,2].  H2  is  a 
traditional  fuel  gas,  however,  high  purity  FI2  is  expensive  and  its 
cost  will  impede  broad  application  of  the  SOFC.  Coal  derived  syngas 
represents  a  more  economical  option  given  the  abundance  of  the 
fuel  as  well  as  the  development  of  gasification  technology.  It  con¬ 
sists  primarily  of  hydrogen  and  CO  with  significant  water  vapour 
and  some  levels  of  C02  and  other  minority  species  [31  ]. 

It  is  commonly  assumed  that  only  H2  oxidation  contributes  to 
electrochemical  power  generation  while  CFI4  is  reformed  to  car¬ 
bon  monoxide,  which  is  then  converted  to  C02  and  H2  through  the 
water-gas  shift  reaction  [3,20,22,32-34].  Consequently,  the  steam 
reforming  reaction  for  methane,  the  water-gas  shift  reaction  and 
electrochemical  reactions  occur  simultaneously  in  the  cell,  and  can 
be  summarized  as  follows: 

CH4  +  H20  -*  CO  +  3H2(reforming)  (1) 

CO  +  H20  -*  C02  +  H2  (water  gas  shift)  (2) 

H2  + 1  /202  -►  FI20  (overall  cell  reaction)  (3) 

Although  a  SOFC  system  may  consist  of  several  repeating  elec¬ 
trochemical  cells  to  reach  technologically  significant  power  levels, 
models  of  SOFCs  are  usually  developed  for  the  smallest  unit-cell 
that  is  assumed  to  represent  the  operating  conditions  of  the  whole 
stack  [23,27].  Thus  a  single  SOFC  cell  simulation  model  is  taken  as 
representative  that  accounts  for  internal  reforming  and  water-gas 
shift  equilibrium,  electrochemical  polarizations  and  the  associated 
heat  generation,  mass  transfer  via  cell  reactions,  and  overall  energy 
balances  [23,26,27].  This  representation  can  be  readily  constructed 
as  quantities  such  as  stack  voltage  and  stack  power  which  are  scaled 
versions  of  single-cell  voltage  and  power. 


3AA.  Mass  balances 

In  terms  of  the  chemical  and  electrochemical  reactions  that 
occur  within  the  anode  and  cathode  of  the  SOFC,  the  first  step  is 
to  calculate  the  equilibrium  gas  composition  at  the  fuel  and  air 
channel  exits. 

The  amount  of  hydrogen  consumed  in  the  fuel  cell  reac¬ 
tions,  hH2-c  (mols-1),  is  related  to  the  current  by  Faraday’s  law 
[22,35,36]: 


where  i  denotes  the  current  density,  A  represents  the  surface  area 
of  the  interconnect  plate  (assuming  the  interconnect  plates  have 
the  same  area),  ne  is  the  number  of  electrons  transferred  in  reac¬ 
tion,  and  F=  96,485  C  mol-1  is  Faraday’s  constant.  For  a  known 
fuel  utilization  factor,  Up  the  amount  of  hydrogen  supplied,  fiH2-s 
(mols-1),  is  given  by 


_  Hh2-c  _  iA 

" h2~ s  _  ~Uf  ~  njvj 


(5) 


The  molar  flow  rate  of  the  fuel  stream,  hf.in  (mol  s-1 )  needed  to 
produce  the  required  amount  of  hydrogen  is  thus 


^  _  hu2-s  _  iA  ^ 

n/'in  ~~  ~  neFUfxfc  1  ] 

where  XfC  is  the  number  of  moles  of  Flydrogen  produced  by  1  mol 
of  fuel  and  it  can  be  calculated  according  to  the  composition  of  the 


fuel  as  XfC  =  Xh2  +  xCo  +  4xCh4-  For  a  known  fuel  gas  composition  Xj, 
its  individual  molar  flow  rate  is: 


iAx  ■ 

n/-i„0)  =  >V-in 


(7) 


where  j  are  the  components  in  the  fuel  stream,  i.e.,  j  =  {H2,  CO,  CFI4, 
C02,H20,N2}. 

An  amount  of  steam  equivalent  to  twice  the  amount  needed  for 
the  reforming  and  water-gas  shift  reactions  is  supplied  in  order  to 
avoid  carbon  deposition  [22].  The  molar  flow  rate  of  steam  needed, 
fiH2o  (mols-1)  is  thus 


Hh2o  =  (n/_in(CO)  +  hf_in(C H4)  X  2)  X  2  (8) 

The  molar  flow  rate  of  additional  steam  supplied,  hADD  (mol  s_1 ) 
is  thus  given  by 


«ADD  =  Hh2o  -  *V-in(Fl20) 

=  (H/-in(CO)  +  hf_in(CH4)  X  2)  X  2  -  hf. in(H20)  (9) 

Therefore  the  total  molar  flow  rate  of  the  fuel  stream  entering  the 
SOFC,  hf  (mol  s-1 ),  becomes 

hf  =  hf_in  +  hADD  (10) 

For  known  conversions  for  the  reforming  and  water  shift  reac¬ 
tions,  the  component  flow  rates  in  the  fuel  exit  stream  are  given 
as 


*V-out(H2)  =  ^h2-s  x  (1  -  Uf)  (11) 

n/-out(CO)  =  0  (12) 

t*/-out(C02)  =  h/-in(CO)  +  h/-_jn(C02)  +  hf_  jn(CH4)  (13) 

n/-out(CH4)  =  0  (14) 

n/-out(H20)  =  ^  +  nH2-s  x  Uf  (15) 

n/-out(N2)  =  n/.in(N2)  (16) 

Given  a  known  inlet  composition,  the  molar  flow  rates  for  the 
air  stream  are: 

Ur/ 2 

*Vin(°2)  =  ^H2-s  x  -j-j—  (17) 

C >a 

na-in(N2)  =  na_in(02)x  57?  (18) 

ha- out(02)  =  ha_in(02)  x  (1  -  Ua)  (19) 

na-out(N2)  =  na_in(N2)  (20) 


where  Ua  is  the  air  utilization  factor,  the  subscript  f  refers  to  the 
fuel  stream,  the  subscript  ‘a’  to  the  air  stream,  the  subscript  ‘in’  to 
the  fuel  cell  inlet  and  the  subscript  ‘out’  the  fuel  cell  outlet. 


3 A. 2.  Electrochemical  descriptions 

For  a  fuel  cell,  its  theoretical  open-circuit  voltage  depends  on 
the  gas  composition  and  temperature  at  the  cathode  and  anode, 
which  can  be  determined  by  the  Nernst  equation  given  as  follows 
[9,23,24,36]: 


-Ag(r,p) 

neF 


(21) 


Note  here  the  molar  Gibbs  free  energy  change  for  the 
fuel  cell  reaction  depends  dramatically  on  the  temperature  (T) 
and  partial  pressures  of  reactants  (p),  i.e.,  A g(T,  p )  =  Ag°(T)  - 
RT  \n(pu2PQ^ /Ph2o)>  where  ft  =  8.314Jmol-1  K-1  is  the  universal 
gas  constant,  Ph2,  Po2>  and  Ph2o  are  the  partial  pressures  of  reac¬ 
tants  H2,  02,  and  H20,  respectively.  Note  that  Ag°(T)  =  Ah0  -  T  As0 
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stands  for  the  molar  Gibbs  free  energy  change  at  p0  =  1  atm  which 
also  depends  on  temperature  [9,23,24,37]. 

While  the  Nernst  potential  represents  the  maximum  theoretical 
reversible  voltage  of  a  fuel  cell  at  the  given  conditions,  the  voltage 
of  an  operating  cell  is  generally  lower  than  this.  As  current  is  drawn 
from  a  fuel  cell,  the  cell  voltage  falls  due  to  internal  resistances  and 
overpotential  losses.  These  losses  are  common  to  all  types  of  fuel 
cells  and  cannot  be  eliminated  [9,24,32].  Here,  three  types  of  polar¬ 
izations,  i.e.,  activation,  ohmic  and  concentration,  are  considered 
and  calculated  through  Eqs.  (22)-(30). 


(1)  Activation  overpotentials  depend  on  the  kinetics  of  the  elec¬ 
trochemical  reactions  occurring  at  the  anode  and  cathode  [2]. 
According  to  the  general  Butler-Volmer  equation  [24,32],  the 
respective  activation  overpotentials  of  the  anode  and  cathode 
can  be  calculated  as  [24,38,39]: 

frsmlr  (sb)  <22> 

W-^inh-’(^)  (23) 

where  i0  a/c  denotes  the  anode/cathode  exchange  current  density 

[32,37,40]. 


Table  1 

Operating  conditions  and  major  parameters  [8,16,20,22,34,37,40]. 


Parameter 

Symbol 

Value 

Ambient  temperature  (I<) 

To 

298 

Operating  pressure  (atm) 

Po 

1 

Fuel  utilization 

Uf 

0.8 

Air  utilization 

Ua 

0.2 

Number  of  electrons 

ne 

2 

Anode  exchange  current  density  (Am-2) 

io,a 

6500 

Cathode  exchange  current  density  (Am-2) 

io,c 

2500 

Limiting  current  density  (Am-2) 

h 

9000 

Anode  thickness  (|jim) 

La 

500 

Anode  conductivity  constants 

CiaJ  6.2a 

95xl06;  -1150 

Cathode  thickness  (|Jim) 

Lc 

50 

Cathode  conductivity  constants 

Cic;  c2c 

42xl06;  -1200 

Electrolyte  thickness  (|Jim) 

Le 

10 

Electrolyte  conductivity  constants 

Cie;  c2e 

3.34  x  104;  -10,300 

Interconnect  thickness  (cm) 

f-int 

0.3 

Interconnect  conductivity  constants 

Clint!  C2int 

9.3  x  106;  -1100 

Compressor  isentropic  efficiency 

Pc 

90% 

Turbine  isentropic  efficiency 

Pe 

93% 

Air  blower  power  consumption  factor 

Pab 

10% 

RT 

^con c,a  —  ^  1 


RT 
neF x 


In  1 


(29) 

(30) 


(2)  Ohmic  overpotentials  are  caused  mostly  by  resistance  to  con¬ 
duction  of  ions  (through  the  electrolyte)  and  electrons  (through 
the  electrodes  and  current  collectors)  and  by  contact  resis¬ 
tance  between  cell  components  [1,2].  The  Ohmic  losses  in  the 
present  study  are  simulated  by  means  of  resistivity,  obtained  by 
the  experimental  temperature-dependent  correlations  [38,39]. 
Such  resistances  are  summated  as  follows  assuming  a  series 
electrical  scheme: 


’'ohm 


=  lYRk  =  =  Y~  =  i(— +  —  +  —  +  — ) 

^  2-^  OkA  Gk  \Oe  Ga  Gc  Gint ) 


(24) 


where  iL,a/c  denotes  the  limiting  current  density  of  the 
anode/cathode. 

Consequently,  the  terminal  voltage  of  the  SOFC  can  be  obtained 
as  follows: 

RT 

V  =  E-Vact-Vohm-Vconc  =  E-—d,  (31) 

Her 

where  d\  =  2  sinh- 1  [i/(2i0>a)]  +  2  sinh-1  [i/(2i0 c)]  -ln(l  -ilh,a)-  In 
(1  -iliL,c)  +  UneFI{RT)]{LelGe  +  LalGa+LclGc  +  LintlGint).  The  electro¬ 
chemical  performance  of  the  fuel  cell  is  described  by  plotting  the 
terminal  voltage  versus  the  current  density  and  operating  temper¬ 
ature  as  shown  in  Fig.  3. 


where  Rk  represents  the  resistance,  Lk  the  thickness,  Ak  the  area, 
and  ak  denotes  the  electronic  conductivity  of  the  anode,  cathode, 
interconnect  and  the  ionic  conductivity  of  the  electrolyte: 


<?e  =  C\e exp  (~j^) 

(25) 

*  =  ^exp(^) 

(26) 

(27) 

<rint  =  %^exp  (^) 

(28) 

where  Cie-C2int  are  constants  defined  in  Refs.  [38,39]  and  shown 
in  Table  1.  The  subscripts  a,  c,  e  and  int  denote  anode,  cathode, 
electrolyte  and  interconnect,  respectively. 


(3 )  Concentration  overpotentials  are  the  voltage  drop  due  to  mass 
transfer  limitations  from  the  gas  phase  into  and  through  the 
electrode  [24].  Previous  papers  based  on  such  models  showed 
that  the  resulting  limiting  current  density  is  only  slightly 
dependent  on  the  SOFC  temperature  [41  ].  Thus,  the  concentra¬ 
tion  overvoltage  calculation  in  the  present  study  is  simplified 
assuming  a  constant  value  for  the  limiting  current  density 
[7,19,32,40],  evaluated  as  follows: 


3.1.3.  Airblower 

Blowers  play  key  roles  in  fuel  cell  designs.  SOFC  systems  require 
blowers  to  provide  motive  force  to  the  incoming  atmospheric  air,  in 
order  to  overcome  the  pressure  drop  in  the  fuel  cell  stack  and  drive 
the  air  through  the  system.  The  electrical  power  required  to  drive 
this  component  is  typically  one  of  the  largest  parasitic  loads  for  the 
fuel  cell  system,  and  one  that,  if  not  carefully  designed  to  meet  the 
external  power  demand,  can  lower  the  overall  system  efficiency 
[22,23,28,42]. 

A  simple  blower  model  is  derived  as  follows  to  determine  the 
power  required  for  blowing  air  in  the  fuel  cell  section: 

Pab  —  Pfc^lab  (32) 

where  rjab,  the  power  consumption  factor,  is  defined  as  a  ratio  of 
the  power  provided  by  the  fuel  cell  itself  for  air  blowing  to  the  total 
amount  of  power  generated  by  the  SOFC,  and  its  value  is  normally 
no  more  than  20%  [42-44].  Obviously  the  lower  value  of  r)ab ,  the 
less  power  supplied  to  the  air  blower  as  parasitic  losses,  and  the 
more  net  power  this  hybrid  system  would  generate,  resulting  in  a 
better  hybrid  performance. 

3.1.4.  Combustor 

Because  only  part  of  the  fuel  can  be  oxidized  in  the  SOFC,  an 
afterburner  is  needed  to  combust  the  residuals  and  produce  addi¬ 
tional  thermal  energy  for  use  elsewhere  in  the  system.  Therefore, 
the  role  of  the  combustor  in  the  hybrid  system  under  study  is  to 
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burn  the  non-reacted  hydrogen  coming  out  on  the  anode  side  of 
the  fuel  cell  with  the  non-reacted  oxygen  exiting  the  cathode  side. 

Using  a  mass  and  energy  balance  on  the  combustor,  a  model 
is  developed  to  determine  the  flow-rate  and  temperature  of  the 
combustor  exit  stream,  given  known  values  for  the  inlet  streams: 


he- out(H2)  —  0 

(33) 

hc-out(CO)  =  0 

(34) 

he- out(C02)  =  h/_out(C02) 

(35) 

hc-out(CH4)  =  0 

(36) 

hc-out(H20)  =  n^_out(H20)  +  hj_  out(^2) 

(37) 

hc-out(N2)  =  h/-out(^2)  +  ha_0ut(N2) 

(38) 

nc-out(02)  =  na.out(02)  hf-outOh) 

(39) 

Given  all  flow-rates,  the  temperature  of  the  combustor  exit  stream 
(Tc)  can  be  further  determined  by  solving  the  energy  balance  equa¬ 
tion  of  the  combustor: 


^/-out(7/-out)  +  ^a-out(^a-out)  —  ^c-out(,Pc)  (40) 

where  T/-_out  =  Ta-0ut  =  T+  AT/2.  Eq.  (40)  implies  that  the  combustor 
exit  temperature  Tc  is  a  function  of  the  fuel  cell  temperature  T. 


3.1.5.  Energy  balances 

Assuming  ideal  gas  behavior,  the  enthalpy  of  a  stream  (W),  H, 
is  calculated  as  a  function  of  the  molar  flow  rates  and  temperature 
[22] by 


H  =  JfnMT)  = 


(41) 


where  the  molar  enthalpy  of  each  species  i  in  the  stream,  frj(T), 
is  given  as  a  function  of  the  local  temperature,  with  A hf  being 
the  standard  enthalpy  change  of  formation  of  species  i  (J  mol-1 ) 
and  cp  i  denoting  the  heat  capacity  of  component  i  (J  mol-1  K-1 ). 
Detailed  polynomial  fitting  functions  for  the  enthalpies,  specific 
heats,  entropies  and  Gibbs  free  energy  values  are  evaluated  based 
on  temperature,  according  to  the  thermodynamic  property  data 
available  from  the  literature  [45-47]. 

Therefore,  the  total  enthalpy  change  for  the  fuel  cell  section  can 
be  determined  as: 

iA 

AH  =  (  —  AHpreheater)  +  AH$oFC  =  — —  (  —  Afapreheater  +  A^SOFc) 
iA 

=  lCFAh  (42) 

where  AHpreheater  denotes  the  enthalpy  change  of  the  endother¬ 
mic  preheating  process,  and  AHSofc  represents  the  enthalpy 
change  for  the  overall  exothermic  electrochemical  reaction. 
Ah  =  -  Ahpreheater +  A/isofc  while  Ahprejlea|-0r  and  AhsoFC  ^re 
expressions  given  in  Appendix  A. 

Considering  the  electrochemical  irreversible  losses  yields  the 
rate  of  the  total  entropy  production  of  the  SOFC  as 


s 


Fig.  2.  The  T-S  diagram  of  an  irreversible  Brayton  cycle. 


where  Pab  accounts  for  the  parasitic  electrical  consumption  in  the 
air  blower  which  is  covered  by  the  electricity  generated  by  the  fuel 
cell.  As  for  the  efficiency,  it  is  defined  as  the  power  output  divided 
by  the  total  energy  input: 


Pfc  =  -A  g-RTd1 
-AH  -Ah(l  +  rjab) 


(45) 


Note  that  the  simulation  of  a  single  cell  unit  is  taken  as  repre¬ 
sentative  of  the  entire  stack  performance. 


3.2.  Gas  turbine  and  heat  exchangers 


3.2 A.  Irreversible  Brayton  cycle 

Gas  turbines  are  power  machines  in  which  the  working  fluid 
is  compressed,  heated,  and  expanded  for  the  purpose  of  generat¬ 
ing  power  [48].  The  basic  thermodynamic  cycle  which  describes 
steady-state  and  steady-flow  energy  characteristics  of  gas  turbine 
components  is  known  as  the  Brayton  cycle  [48].  An  irreversible 
Brayton  cycle  is  depicted  in  Fig.  2,  where  process  1-2S  is  an 
isentropic  compression  representing  the  power  used  by  an  ideal 
compressor,  while  process  1-2  takes  into  account  the  non- 
isentropic  nature  of  a  real  compressor.  Process  2-3  is  the  isobaric 
heat  addition  and  process  4-1  is  an  isobaric  heat  rejection.  The 
process  3-4S  is  an  isentropic  expansion  representing  the  power 
output  from  an  ideal  turbine,  while  process  3-4  takes  into  account 
the  non-isentropic  nature  of  a  real  turbine.  Tj(i  =  l,  2,  3,  4)  is  the 
temperature  of  the  working  fluid  in  the  gas  turbine  at  state  points 
1,  2,  3,  and  4.  and  Qi  denote  the  heat  flow  transferred  via  HEX1 
and  HEX2,  respectively.  In  the  present  paper,  a  practical  irreversible 
cycle  consisting  of  states  1-2-3-4-1  is  considered. 

To  describe  the  irreversibility  of  a  practical  gas  turbine  cycle,  the 
isentropic  efficiencies  of  the  compressor  and  turbine  are  introduced 
as  controllable  parameters  [49,50]: 


Be  = 


T2S  -  P\ 
t2-t i 


(46) 


Stot  — 


f(^act  +  ^ohm  +  ^conc) 

T0 


T3-T4 

(43)  Be  =  j3_Tas 


(47) 


where  To  is  the  ambient  temperature.  Based  on  this,  the  net  power 
output  of  the  fuel  cell  can  be  deduced  as  a  function  of  current 
density,  temperature,  partial  pressures,  chemical  composition,  and 
geometric/material  characteristics  as: 


When  T]c  =  1  and  r\e  =  1 ,  P2  =  T2S,  T4  =  T4S,  and  the  irreversible  com¬ 
pression  and  expansion  processes  become  reversible.  According  to 
the  adiabatic  equations  of  an  ideal  gas  and  the  characteristics  of  the 
Brayton  cycle: 


Pfc 


-AG  -  ToStot  -  Pab 


iA(-Ag-RTd 1) 
neT(l  +  T]ai,) 


(44) 


P3  _  T4S 
T2s  Pi 


(48) 
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Fig.  3.  The  fuel  cell  terminal  voltage  varying  with  current  density  (i)  and  operating 
temperature  (T),  where  the  relevant  parameters  are  listed  in  Tables  1  and  2. 


the  turbine  outlet  temperature  can  be  readily  deduced  from  Eqs. 
(46)-(48)  as: 


T4  —  (1  —  rje)T 3  + 


xr)eT 1 _ 

ric  +  x{\  -77c)(7i/r3) 


(49) 


where  x  =  T3/T2  is  defined  as  the  temperature  ratio  of  the  isobaric 
process  2-3,  namely,  the  ratio  of  turbine  inlet  temperature  to  com¬ 
pressor  exit  temperature.  Note  that  instead  of  using  turbine  inlet 
temperature  as  an  operating  parameter  [8,21,51],  this  study  has 
introduced  x  as  a  major  controllable  variable  to  determine  the 
contribution  of  GT  cycle  components  on  the  entire  system  perfor¬ 
mance. 


3.2.2.  Heat  exchangers 

As  illustrated  in  Fig.  1,  three  heat  exchangers  are  utilized  in  the 
system,  all  of  which  are  modelled  as  flat  plate  counter  flow  heat 
exchangers  as  presented  in  Refs.  [9,52].  The  primary  heat  exchanger 


(HEX1),  which  plays  a  vital  role  in  the  hybrid  cycle,  is  used  to 
recover  the  heat  of  the  high  temperature  fuel  cell  exhaust  to  drive 
the  bottoming  gas  turbine.  Moreover,  fuel  and  air  pre-heaters  are 
required,  as  the  fuel  cell  does  not  tolerate  a  gas  supply  at  low  tem¬ 
peratures  because  of  the  excessive  cooling  and  consequent  thermal 
stresses  that  cold  streams  would  cause.  Therefore,  two  additional 
heat  exchangers  (HEX2  and  HEX3)  are  further  adopted  to  fully  recy¬ 
cle  the  waste  heat  of  the  system  and  preheat  the  fuel  cell  incoming 
reactants  in  two  steps.  The  heat  transfer  processes  involved  in  the 
two  pre-heaters  are  assumed  to  be  ideal  for  the  sake  of  simplicity. 

As  a  main  source  of  irreversibility,  the  loss  caused  by  the 
finite-rate  heat  transfer  between  the  SOFC  and  GT  subsystems  is 
considered.  Once  the  finite  heat  transfer  rate  is  taken  into  account, 
the  performance  of  the  gas  turbine  is  closely  dependent  on  heat- 
transfer  laws.  According  to  Newtonian  heat-transfer  [53]  and  the 
expression  of  log  mean  temperature  difference  (LMTD)  [54],  the 
heat  transfer  rates  for  the  counter-flow  heat  exchanger,  and  Qi , 
can  be  expressed  as 

=  (50, 


Qi  =  mCp{T4  -  Ti)  = 


U2A2^+T0-^-rle)T3-vT,] 
ln{l /[( 1  -  r]e)T3  +  vT\  -  Tq]} 


(51) 


where  u  =  (Tc- T3)/(T-  AT/2  -  I3/x),  v  =  xr\el\r\c  +  x(l  - 

^c)(Ti/T3)],  m  and  cp  denote,  respectively,  the  mass  flow  rate 
and  heat  capacity  at  constant  pressure  of  the  working  fluid  in  the 
gas  turbine,  Lfl  and  U2  are  the  heat  transfer  coefficients  of  HEX1 
and  HEX2,  A\  and  A2  represent  the  corresponding  heat  transfer 
surface  areas. 

Apart  from  the  irreversibility  existing  in  heat  exchangers,  the 
other  major  loss  is  due  to  heat  lost  from  the  SOFC  directly  to  the 
environment.  Heat  loss  is  modelled  based  on  the  temperature  dif¬ 
ference  between  the  fuel  cell  and  ambient  conditions  as  follows 
[7,19,55]: 


Qloss  =  XA/(T  —  Tq) 


(52) 


Fig.  4.  The  optimized  GT  temperatures  varying  with  the  operating  temperature  (T)  and  temperature  difference  (AT)  of  the  fuel  cell,  where  the  parameters  i  =  2000  Am-2, 
x  =  1.5,  mi  =1.1,  m2  =  -1.2  x  10-7,  and  m3  =0.01  are  chosen.  The  values  of  other  parameters  are  listed  in  Tables  1  and  2. 
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where  I<  represents  the  effective  convective  and/or  conductive 
heat-leak  coefficient,  and  A/  denotes  the  effective  heat-transfer 
area.  Note  that  this  equation  captures  the  effects  of  fuel  cell  temper¬ 
ature  changes.  Based  on  this,  the  heat  flow  rates  transferred  from 
the  SOFC  to  the  gas  turbine  can  be  obtained: 

Qh  =  -  AH  -  Pfc  -  Qoss  =  - AH(1  -  rift)  -  KAi(T  -  T0)  (53) 

Therefore  the  efficiency  of  the  gas  turbine  cycle  can  be  expressed 
as 

.  Qt  ,  mcp(T4-TQ  (1  —  ^c)T3 +FT] -Tt 
(j,  mcp(T3-T2)  T3-T3/x 

(54a) 


or 


entire  SOFC-GT  hybrid  system  can  be  obtained: 


*7 hybrid  — 


^hybrid 

-AH 


l Ifc  + 


% 

-AH 


—  rife  + 


Qh 

-AH 


•  dgt 


=  rife  +  rjgt 


1  -Vfc  - 


m3(T  -  Tp)- 
ini2  A  h 


(59) 


^hybrid  ~  AH  •  ^hybrid 


iA  Ah 
neF 


rife  +  rjgt 


1  -rife  ~ 


m3 (T  -  T0 


irri2  A  h 


n 


(60) 


It  is  clearly  seen  from  Eqs.  (59)  and  (60)  that  the  efficiency 
and  power  output  of  the  hybrid  system  are  closely  dependent 


„  ,  Qi  ,  AIi/1(w-r1xi  +  r0-w)[rc-r  +  Ar/2-r3(i-i/x)] 

lint  —  1  —  — —  —  1  —  — - 7 - 

Qfc  Qh  {(Ui/l/2)(w  -  T,  )[TC  -  T  +  AT/2  -  r3(  1  -  l/x)]ln[l/(w  -  T0)]  +  T3(l  -  l/x)(l  +  T0  -  w)lnu} 

where  w  =  (1  -  rje)T3  +vT^,  and  =A\  +A2  is  the  total  heat  trans¬ 
fer  surface  area.  Comparison  of  Eq.  (54a)  with  Eq.  (54b)  yields  the 
following  relation: 

_ AhH1T3(l  -  l/x)(l  +  Tp  -  w)[Tc  -  T  +  AT/2  -  I3(l  -  1/x)] _ =  1 

QftKL/i/^Xw  -  7i )[TC  -  T  +  AT/2  -  T3(l  -  l/x)]ln[l/(w  -  T0)]  +  T3(l  -  l/x)(l  +  T0  -  w)lnu} 


(54b) 


(55) 


Eq.  (55)  indicates  that  Ti  is  a  function  of  T3.  Using  the  solution  of  Ti 
in  Eq.  (55),  we  can  eliminate  T\  from  Eqs.  (54a)  and  (54b).  There¬ 
fore,  the  gas  turbine  efficiency  becomes  only  a  function  of  T3 ,  Ah  and 
Qh  with  other  given  parameters,  i.e.,  rjgt  =  rjgt{T3 ,  Ah,  Qh).  Accord¬ 
ingly,  when  Qh  (the  amount  of  heat  released  from  the  fuel  cell  to 
the  gas  turbine)  and  Ah  (the  total  heat  transfer  area  of  the  gas  tur¬ 
bine  cycle)  are  given,  using  Eqs.  (54a)  and  (54b)  and  their  extremal 
conditions  (drjgt  /  dT3)A^  ^  =  0,  it  can  be  proved  that  when  the  gas 
turbine  efficiency  attains  its  maximum,  the  relationship  between 
the  turbine  inlet  temperature  T3  and  compressor  inlet  temperature 
T\  is  determined  by  the  following  equation: 


on  the  parameters  related  to  the  various  irreversible  losses,  which 
include  the  irreversibilities  within  the  fuel  cell  itself  and  orig¬ 
inating  from  heat  transfer  due  to  convection/conduction  in  the 
combined  SOFC-GT  power  plant.  The  model  is  based  on  the  physics, 
chemistry  and  electrochemistry  that  govern  the  entire  system, 
and  is  developed  by  resolving  conservation  of  mass  and  energy 
principles. 


[(1  -  l/x)(7i  -D3)DS+D4D5D3+D4D3{TX  -D3)]D6+DaD5(T,  -D3)(D13 +D14)  _  Ti 

[D4D5(D9  +  1)  +  D4D9(Tt  -  D3)]D6  -  D4D5D15(T1  —  d3)  t3 


(56) 


where  D\  -D\3  are  expressions  given  in  Appendix  A. 

Substituting  the  solution  of  T[  from  Eq.  (56)  into  Eq.  (55)  enables 
us  to  determine  the  relationship  of  the  optimized  T3  varying  with 
T  as  illustrated  in  Fig.  4c,  which  shows  that  the  optimized  value  of 
turbine  inlet  temperature  is  a  monotonically  increasing  function  of 
the  fuel  cell  operating  temperature  for  other  given  parameters. 

Following  the  above  procedure,  the  optimized  GT  efficiency  rjgt 
becomes  only  a  function  of  the  fuel  cell  temperature  T  for  a  set 
of  given  parameters.  Combining  the  above  equations  yields  the 
optimum  efficiency  and  power  output  of  the  gas  turbine  as: 


4.  Parametric  optimization  analysis 

In  order  to  determine  the  sensitivity  of  the  overall  hybrid  per¬ 
formance  to  the  typical  component  parameters  and  to  explore 
the  optimum  system  operation,  numerical  simulations  of  this 
combined  SOFC-GT  power  plant  are  carried  out  using  MAT- 
LAB.  The  simulation  model  is  based  on  a  number  of  parameters, 
dramatically  affecting  the  results  coming  out  from  the  code. 


(w  -  Ti  )(1  +  7b  -  w)[Tc  -  T  +  AT/2  -  T3(l  -  1/x)] 

\im2  Ah(l  -  rjfc)  -  m3(J  -  T0)]{mi(w  -  T^Tc  -  T  +  AT/2  -  T3(l  -  l/x)]ln[l/(w  -  T0)]  +  T3(l  -  l/x)(l  +  T0  -  w)lnu} 


(57) 


Pgt  —  Qh^lgt  —  rjgt 


-  iA  Ah 
neF 


( 1  rjfc )  + 


Am3 

neFiri2 


(T  -  Tq) 


(58) 


where  m\  =  Ui/U2,  m2  =  -A/(neFAhU^)  and  m3  =/<A//(UiA/l).  Note 
here  temperatures  T3  and  are  determined  by  simultaneously 
solving  Eqs.  (55)  and  (56).  Eqs.  (57)  and  (58)  indicate  that  the  opti¬ 
mized  efficiency  and  power  output  of  the  SOFC-driven  gas  turbine 
are  closely  dependent  on  the  design  and  operating  parameters  of 
the  upstream  fuel  cell  subsystem. 


3.3.  Ffybrid  efficiency  and  power  output 

Based  on  the  joint  solution  of  Eqs.  (53)-(58),  the  following 
expressions  of  the  optimized  efficiency  and  power  output  for  the 


The  model  easily  allows  changes  and  parameterization  of  any  vari¬ 
able  in  it.  For  example,  the  fuel  cell  current  density  (z),  operating 
temperature  (T),  temperature  difference  (AT),  fuel  utilization  factor 
(Uf),  air  blower  power  consumption  factor  ( rjab ),  the  GT  temper¬ 
ature  ratio  (x)  and  isentropic  efficiencies  {r]c  and  rje ),  as  well  as 
the  synthesized  parameters  related  to  the  finite-rate  heat  transfer 
via  heat  exchangers  (mi  and  m2\  etc.,  are  treated  as  independent 
design  parameters  that  can  individually  be  varied.  For  purpose 
of  validation,  the  overall  performance  is  predicted  by  the  present 
model  based  on  the  validated  data  adopted  from  previous  literature 
[8,20,24,26,38,41,44,56].  Values  of  those  parameters  are  summa¬ 
rized  in  Tables  1  and  2,  and  are  used  as  constants  throughout  the 
study  unless  mentioned  specifically. 


Table  2 

Coal  syngas  composition  [52]. 
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Components 

Mole  fractions 

h2 

0.293 

CO 

0.287 

co2 

0.118 

ch4 

0.000 

n2 

0.030 

h2o 

0.272 

In  particular,  some  operating  constraints  have  been  considered 
in  order  to  determine  a  more  generalized  system.  For  example,  the 
fuel  cell  temperature  rise,  i.e.,  the  temperature  difference  between 
the  SOFC  inlet  and  outlet,  must  remain  below  200  K  to  minimize 
thermal  stresses  in  a  practical  fuel  cell  since  its  tri-layer  materials 
are  ceramics  which  are  not  capable  of  withstanding  the  thermal 
stresses  associated  with  larger  temperature  gradients  [23].  Thus, 
a  temperature  difference  greater  than  200  K  is  excluded  from  the 
practical  design  range.  As  for  the  range  of  fuel  utilizations  con¬ 
sidered,  although  fuel  cell  efficiency  continuously  increases  with 
increasing  fuel  utilization,  utilizations  of  90%  and  above  are  not 
practically  achievable  due  to  imperfect  fuel  distribution  and  the 
need  for  a  margin  of  safety  [23].  Therefore  the  average  fuel  uti¬ 
lization  selected  for  most  of  the  simulation  in  this  paper  is  80%. 
Moreover,  as  a  synthesized  parameter  to  measure  the  effect  of  heat 
loss  irreversibility  on  the  system  performance,  different  values  of 
m3  correspond  to  different  quantities  of  energy  released  by  the  sys¬ 
tem  as  heat  loss  to  the  environment.  A  value  of  0.01  for  m3  is  found 
to  result  in  approximately  7-9%  heat  loss  with  respect  to  the  total 
energy  input  of  the  system  at  typical  operating  conditions,  which 
is  quite  reasonable  according  to  previous  studies  [23].  Hence  this 
value  is  chosen  for  m3  throughout  the  simulation  analysis. 

Using  Eqs.  (59)  and  (60),  as  well  as  the  solutions  of  T3  and 
Ti  which  are  determined  by  the  simultaneous  Eqs.  (55)  and  (56), 
enables  us  to  plot  the  efficiency  and  power  density  of  the  opti¬ 
mized  SOFC-GT  system  varying  with  the  fuel  cell  current  density 
and  operating  temperature,  as  presented  in  Fig.  5,  where  the  power 
density  hybrid  =  ^hybrid /A  The  meshes  are  generated  over  a  spec¬ 
ified  temperature  domain  of  873-1273  K  and  a  current  density 
domain  of  1800-2200 Am-2  because  the  SOFC  generally  oper¬ 
ates  within  these  ranges  [8,51].  It  is  interesting  to  observe  that 
the  hybrid  efficiency  passes  through  a  maximum  value  when 
the  SOFC  temperature  is  varied,  while  the  hybrid  power  density 
increases  monotonically  with  increasing  temperature.  Therefore 
there  exists  an  optimum  T0pt  which  denotes  the  corresponding 
temperature  when  the  hybrid  efficiency  attains  its  maximum  r]max. 
This  is  different  from  the  cases  reported  in  Refs.  [21,23]  that  indi¬ 
cated  monotonically  increasing  efficiency  of  the  hybrid  system  with 
increasing  fuel  cell  temperature.  On  the  other  hand,  when  the  cur¬ 
rent  density  of  the  fuel  cell  is  increased  from  1800  to  2200  Am-2, 
the  hybrid  efficiency  will  decrease  monotonically  while  there  exists 
a  maximum  power  output  Pmax  with  iopt  being  the  corresponding 
optimum  current  density.  Note  here  the  behavior  of  system  effi¬ 
ciency  is  consistent  with  most  of  the  cases  for  SOFC-GT  hybrids 
[23,57-59].  Furthermore,  Topt  will  decrease  with  the  increase  in  fuel 
cell  current  density,  while  iopt  increases  as  the  fuel  cell  temperature 
is  increased. 

The  advantage  of  the  hybrid  configuration  can  be  better  illus¬ 
trated  by  comparing  its  performance  with  that  of  the  SOFC  only  or 
GT  only  system.  As  shown  in  Fig.  6,  under  typical  operating  con¬ 
ditions  the  maximum  hybrid  power  density  is  2300  Wm-2,  about 
68%  of  that  is  from  the  SOFC  and  32%  is  from  GT.  The  maximum  effi¬ 
ciency  of  the  hybrid  system  is  57%,  which  is  15%  higher  than  that 
of  the  SOFC  only  system  and  19%  higher  than  the  GT  only  system. 


Fig.  5.  The  system  efficiency  (a)  and  power  density  (b)  of  the  hybrid  SOFC-GT 
varying  with  the  operating  temperature  (T)  and  current  density  (i)  of  the  fuel  cell, 
where  the  parameters  AT=  100 K,  x  =  1.5,  m\  =  1.1,  m2  =  -1.2  x  10-7,  and  m3  =0.01 
are  chosen.  The  values  of  other  parameters  are  listed  in  Tables  1  and  2. 


Accordingly,  having  established  the  baseline  performance,  more 
selected  system  parameters  will  be  varied  and  the  sensitivity  of 
these  parameters  to  the  hybrid  performance  will  be  evaluated  in 
the  following  sections. 

4.  I .  Influence  of  temperature 

In  a  hybrid  system,  temperatures  of  different  components  are 
strongly  coupled.  Therefore,  knowledge  of  the  prime  temperature- 
related  parameters  is  required.  There  are  four  major  temperature 
parameters:  the  turbine  inlet  temperature  (I3),  the  fuel  cell  oper¬ 
ating  temperature  (T),  the  GT  temperature  ratio  (x),  and  the 
temperature  difference  across  the  fuel  cell  (AT).  Only  three  of 
them  can  be  selected  as  independent  parameters  and  the  remain¬ 
ing  parameter  is  obtained  from  the  analysis  [21].  In  the  present 
study,  T,  AT  and  x  are  set  as  the  independent  design  parameters 
and  the  turbine  inlet  temperature,  i.e.,  T3,  is  determined  from  the 
optimization  procedure  as  described  above. 

To  illustrate  how  the  SOFC  temperature  affects  the  GT  cycle, 
Fig.  4  has  been  plotted  to  show  the  sensitivity  of  the  GT  temper¬ 
atures  to  the  operating  temperature  and  temperature  difference 
of  the  fuel  cell.  As  shown  in  Fig.  4c,  increasing  the  fuel  cell  tem¬ 
perature  (T)  tends  to  increase  the  turbine  inlet  temperature  (T3). 
For  a  given  GT  temperature  ratio  (x),  this  would  of  course  result  in 
an  increased  compressor  exit  temperature  (T2)  as  seen  from  Fig.  4b. 
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T(K)  T(K) 

Fig.  6.  The  efficiency  and  power  density  comparisons  between  SOFC,  GT  and  hybrid  system,  where  the  parameters  i  =  2000  A  m- 2 ,  AT=  1 00  K,  x  =  1 .5,  m\  =  1 . 1 ,  m2  =  - 1 .2  x  1 0-7 , 
and  m3  =0.01  are  chosen.  The  values  of  other  relevant  parameters  are  listed  in  Tables  1  and  2. 


Meanwhile,  it  is  found  that  both  the  turbine  outlet  temperature  (T4) 
and  the  compressor  inlet  temperature  (7^ )  will  first  decrease  and 
then  increase  with  an  increasing  T,  which  can  be  further  explained 
by  Eqs.  (49),  (55)  and  (56).  On  the  other  hand,  the  variation  in  the 
fuel  cell  temperature  difference  seems  to  have  a  relatively  small 
effect  on  the  GT  temperatures.  As  shown  in  Fig.  4,  both  T2  and  T3  will 
increase  with  an  increasing  AT,  while  both  T\  and  T4  will  decrease 
as  AT  is  increased,  although  these  variations  are  not  very  obvious. 

The  effect  of  SOFC  operating  temperature  on  the  hybrid  per¬ 
formance  can  be  further  explained  by  the  curves  in  Figs.  3-6.  As 
previously  shown  in  Fig.  5,  the  hybrid  efficiency  will  pass  through  a 
maximum  value  while  the  power  density  only  increases  monoton- 
ically  and  quickly  with  an  increase  in  T.  This  is  due  to  many  reasons. 
Firstly,  the  fuel  cell  terminal  voltage  with  respect  to  temperature 
rise  has  a  peak  instead  of  a  monotonically  changing  trend  according 
to  Fig.  3.  Hence  for  a  given  current  density,  maximum  values  of  the 
power  density  and  cell  efficiency  are  achievable,  shown  by  the  blue 
curves  in  Fig.  6.  At  the  same  time,  increasing  the  fuel  cell  tempera¬ 
ture  leads  to  a  greater  GT  power  and  efficiency  as  described  by  the 
red  curves  in  Fig.  6.  This  is  mainly  because  the  turbine  inlet  tem¬ 
perature  is  directly  increased  with  the  fuel  cell  temperature  rise, 
as  can  be  readily  seen  from  Fig.  4.  When  the  fuel  cell  operates  at  a 
higher  temperature  its  power  contribution  decreases,  however,  this 
reduction  is  compensated  by  a  larger  amount  of  power  produced 
by  the  GT  at  higher  temperature,  thereby  leading  to  an  increase  in 
the  overall  hybrid  power  output  as  shown  by  the  black  curve  in 
Fig.  6.  In  addition,  increasing  the  fuel  cell  temperature  also  means 
an  increase  in  the  total  energy  input  of  the  hybrid  system,  which 
conclusively  results  in  an  optimal  hybrid  efficiency  as  described  in 
Fig.  5.  (For  interpretation  of  the  references  to  colour  in  this  text,  the 
reader  is  referred  to  the  web  version  of  the  article.) 

On  the  gas  turbine  side,  its  temperature  ratio  (x)  is  another  inde¬ 
pendent  design  parameter.  Concerning  the  sensitivity  of  the  hybrid 
performance  to  the  GT  temperature  ratio  and  fuel  cell  tempera¬ 
ture  difference,  Fig.  7  provides  a  good  example  for  the  comparison 
of  these  two  different  parametric  effects.  According  to  the  sim¬ 


ulation  results,  the  hybrid  efficiency  increases  dramatically  with 
increasing  x,  which  is  mainly  due  to  the  higher  turbine  inlet  temper¬ 
ature  and  subsequently  an  increased  GT  power  contribution.  When 
more  GT  power  is  generated  the  total  power  output  of  the  system 
will  increase  accordingly,  and  consequently  the  hybrid  efficiency 
becomes  higher.  Compared  to  the  effect  of  x,  the  hybrid  efficiency 
remains  almost  constant  when  AT  is  varied.  This  is  quite  different 
from  the  conclusions  of  Refs.  [21,23]  which  indicate  a  huge  effect 
of  AT  on  the  system  performance.  Nevertheless,  the  slight  reduc¬ 
tion  in  the  hybrid  efficiency  with  decreasing  AT  is  consistent  with 
previous  results  [21].  The  main  cause  is  due  to  the  lower  turbine 
inlet  temperature  for  a  smaller  AT  and  the  resulting  reduced  power 
contribution  of  the  gas  turbine. 


Fig.  7.  The  system  efficiency  of  the  hybrid  SOFC-GT  varying  with  the  GT  tempera¬ 
ture  ratio  (x)  and  the  fuel  cell  temperature  difference  (AT),  where  the  parameters 
z'  =  2000 Am-2,  T=  1073 K,  mi  =1.1,  m2  =  -1.2  x  10-7,  and  m3  =  0.01  are  chosen.  The 
values  of  other  relevant  parameters  are  listed  in  Tables  1  and  2. 
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Fig.  8.  The  system  efficiency  of  the  hybrid  SOFC-GT  varying  with  the  compres¬ 
sor  isentropic  efficiency  (rjc)  and  the  turbine  isentropic  efficiency  ( rje ),  where  the 
parameters  i  =  2000  Anr2,  T=  1073 K,  A T=  100K,  x  =  1.5,  mi  =  1.1,  m2  =  -1.2  x  10-7, 
and  m3  =0.01  are  chosen.  The  values  of  other  relevant  parameters  are  listed  in 
Tables  1  and  2. 


4.2.  Influence  ofGT  isentropic  efficiencies 

As  seen  in  Fig.  8,  the  irreversibility  in  the  adiabatic  processes  of 
the  gas  turbine  cycle  has  great  effect  on  the  system  performance. 
The  hybrid  efficiency  will  increase  rapidly  when  either  of  the  isen¬ 
tropic  efficiencies  of  the  compressor  and  turbine  are  increased.  This 
is  mainly  due  to  the  higher  turbine  inlet  temperature  and  subse¬ 
quently  a  raised  GT  power  contribution  that  is  caused  by  an  increase 
in  turbine  and  compressor  efficiencies.  However,  it  shows  that 
the  hybrid  efficiency  increases  more  significantly  with  the  turbine 
isentropic  efficiency  (r]e)  than  the  compressor  isentropic  efficiency 
According  to  the  simulation  results,  hybrid  efficiency  tends 
to  increase  by  around  17%  when  rje  increases  from  0.8  to  1,  and 
this  efficiency  growth  is  9%  higher  than  that  with  an  increasing 
Pc.  It  is  noteworthy  that  such  a  high  sensitivity  of  the  system  per¬ 
formance  to  changes  in  the  GT  expansion  process  has  not  been 
reported  before.  A  better  turbine  design  is  obviously  quite  essential 
to  enhance  the  entire  hybrid  SOFC-GT  performance. 

4.3.  Influence  of  finite-rate  heat  transfer 

The  effect  of  irreversibilities  due  to  finite-rate  heat  transfer  on 
the  system  performance  can  be  best  explained  by  two  parameters, 
m\  and  m2. 

As  the  ratio  of  heat  transfer  coefficients  between  two  heat 
exchangers  (HEX1  and  HEX2),  m\  =  U\IU2  is  only  to  measure  the 
performance  of  the  GT  cycle  and  is  independent  of  the  upstream 
SOFC  subsystem.  Fig.  9  shows  clearly  that  the  predicted  hybrid  effi¬ 
ciency  will  decrease  as  parameter  m\  is  increased.  This  is  because 
an  increasing  m\  will  lead  to  a  lower  turbine  inlet  temperature, 
which  in  turn  results  in  a  decreased  GT  power  contribution.  Thus 
due  to  the  relatively  smaller  contribution  of  the  gas  turbine  power 
the  total  system  power  is  decreased,  and  accordingly  the  hybrid 
efficiency  decreases. 

On  the  contrary,  the  hybrid  efficiency  tends  to  increase  rapidly 
with  increasing  m2 ,  which  is  a  synthetic  parameter  to  measure  the 
irreversibility  of  finite-rate  heat  transfer  in  the  gas  turbine,  as  also 
illustrated  in  Fig.  9.  This  can  be  attributed  to  the  higher  turbine  inlet 
temperature,  and  thus  a  bigger  GT  power  contribution,  resulting 
from  an  increased  m2. 

In  conclusion,  a  higher  value  of  m2  and  a  lower  value  of  m\ 
would  lower  heat  transfer  irreversibilities  and  thus  improve  hybrid 
performance. 
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Fig.  9.  The  system  efficiency  of  the  hybrid  SOFC-GT  varying  with  mi  and  m2,  where 
the  parameters  1  =  2000 Am-2,  T=  1073  K,  AT=  100K,  x  =  1.5,  and  m3  =0.01  are  cho¬ 
sen.  The  values  of  other  relevant  parameters  are  listed  in  Tables  1  and  2. 


4.4.  Influence  of  fuel  utilization  and  electrical  parasitic  losses 

As  depicted  visibly  in  Fig.  10,  fuel  utilization  and  electrical  par¬ 
asitic  consumption  both  affect  the  performance  of  the  combined 
SOFC-GT  system.  First,  the  hybrid  efficiency  tends  to  increase  by 
only  1.5%  when  the  fuel  utilization  factor  (Uf)  increases  from  0.7  to 
0.9.  The  slight  increase  is  due  to  many  reasons.  First  of  all,  when 
operating  at  higher  fuel  utilization  factor,  due  to  decreased  inlet 
fuel  flow  rate,  the  SOFC  works  at  lower  voltage  which  in  turn 
means  reduced  SOFC  power  contribution.  However,  higher  fuel  uti¬ 
lization  would  result  in  higher  turbine  inlet  temperature,  leading 
to  increased  GT  power  contribution.  Hence,  an  enhanced  overall 
hybrid  efficiency  is  expected,  illustrated  in  Fig.  10. 

On  the  other  hand,  an  increase  by  9%  in  the  hybrid  efficiency 
is  predicted  when  the  power  consumption  factor  of  the  air  blower 
(pab)  decreases  from  20%  to  0%,  due  to  the  decreased  parasitic  losses 
in  the  system. 

Comparing  the  influence  of  these  two  factors,  the  magnitude 
for  that  of  the  parasitic  losses  is  much  more  significant.  Therefore, 
more  effort  should  be  put  into  the  development  of  more  efficient 
blowers  to  reduce  the  system  parasitic  losses. 


Fig.  10.  The  system  efficiency  of  the  hybrid  SOFC-GT  varying  with  fuel  uti¬ 
lization  factor  Uf  and  air  blower  power  consumption  factor  r]ab,  where  the 
parameters  i  =  2000  A  nr2,  T=  1073  K,  AT=  100K,  x  =  1.5,  mi  =  1.1,  m2  =  -1.2  x  10“7, 
and  m3  =  0.01  are  chosen.  The  values  of  other  relevant  parameters  are  listed  in 
Tables  1  and  2. 
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5.  Conclusions 

In  this  article,  an  integrated  atmospheric  syngas  fueled  SOFC-GT 
power  plant  has  been  examined.  An  optimization  strategy  apply¬ 
ing  a  generalized  theoretical  model  of  this  hybrid  has  been  carried 
out  through  thermodynamic  system  modelling  to  determine  the 
optimum  integration  between  the  fuel  cell  and  gas  turbine.  The 
governing  equations  of  energy  conversion  for  each  system  compo¬ 
nent  as  well  as  the  entire  system  have  been  developed.  The  overall 
system  performance  and  the  effect  of  varying  design  and  control 
parameters  have  been  predicted  through  a  detailed  parametric 
analysis,  where  the  current  density,  operating  temperature,  fuel 
utilization  factor  and  temperature  gradient  of  the  fuel  cell,  the  isen- 
tropic  efficiencies  and  temperature  ratio  of  the  gas  turbine  cycle, 
as  well  as  the  parameters  related  to  the  finite-rate  heat  transfer 
between  subsystems  are  all  set  to  be  controllable  variables. 

The  theoretical  analyses  have  indicated  that  by  selecting  an  opti¬ 
mal  turbine  inlet  temperature  according  to  the  joint  solution  of  Eqs. 
(55)  and  (56),  the  highest  GT  efficiency  and  the  maximum  power 
output  will  be  achieved,  which  in  turn  implies  an  optimum  inte¬ 
gration  between  the  fuel  cell  and  the  downstream  gas  turbine. 
In  addition,  the  parametric  analyses  have  identified  a  number  of 
essential  performance  characteristics.  For  example,  the  impact  of 
integrating  a  GT  bottoming  cycle  is  found  to  be  significant,  adding 
as  much  as  19%  of  efficiency  to  the  system  at  typical  operating  con¬ 
ditions.  Both  the  hybrid  efficiency  and  power  density  will  change 
dramatically  with  variations  in  the  fuel  cell  current  density  and 
operating  temperature  over  the  ranges  investigated.  Particularly, 
the  simulation  results  revealed  that  an  increase  in  either  compres¬ 
sor  isentropic  efficiency  or  turbine  isentropic  efficiency  will  lead  to 
a  significant  improvement  in  the  overall  efficiency  of  the  system, 


shown  to  give  comparable  performance  predictions.  The  designed 
optimization  strategy  presented  here  provides  valuable  insights 
regarding  design  tradeoffs  and  optimal  integration  considera¬ 
tions  for  hybrid  fuel  cell  power  plants.  To  make  the  performance 
improvements  potentially  achievable,  the  analysis  reveals  that  it 
is  desirable  to  design  such  hybrid  systems  with  the  largest  possi¬ 
ble  GT  temperature  ratio,  the  lowest  possible  parasitic  losses,  and 
moderate  fuel  cell  current  density  and  operating  temperature. 

Finally,  this  analysis  will  be  the  basis  for  a  future  determina¬ 
tion  of  the  system  performance  capability  and  robustness  to  fuel 
concentration  variations  by  considering  more  complex  models.  To 
explore  better  optimization  methodologies,  the  practical  feasibility 
of  such  optimization  procedures  still  needs  to  be  investigated,  and  a 
more  in-depth  analysis  of  the  system  response  to  a  larger  number  of 
decision  variables  than  actually  considered  in  the  present  study  will 
be  part  of  our  ongoing  research.  The  model  would  also  benefit  from 
validation  against  experimental  data  based  around  the  cycles  anal¬ 
ysed  here,  though  such  data  is  lacking  in  the  literature  at  present. 
Note  that  the  performance  characteristics  presented  herein  may 
change  significantly  with  various  system  configurations,  and  this  is 
also  a  topic  which  should  be  considered  for  future  studies. 
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Appendix  A. 

1.  The  expressions  of  Aftpreheater  and  AftSoFC  in  Eq.  (42)  are  given 

as  follows: 


Afrpreheater  =  [xH2ftH2(’E)  +  xCoftco(E)  +  XcH4ftcH4(E)  +  XCo2ftc02(’E)  +  xN2ftN2(E)  +  (2xCo  +  4xCh4  )ftH2o(E)J  /(Lf/Xfc) 
ho2(T)  +  3.7619hN2(r)]  /(2Ua) 

*H2/!H2(r0)  +  XcohcoiTo)  +  XcH4hcH4{T0)  +  Xco2/lc02(ro)  +  *N2/lN2(7o)  +  (2Xco  +  4xCH4)/lH2o(To)]  /( UfXfc ) 

ho2(r0)  +  3.7619fiN2(r0)]  /(2 Ua) 


AhsoFC  = 


1(1  -  Uf)Xfc  ■  Iih2(T)  +  (xco  +Xco2  +Xch4)  '  hco2(J)  +  (xco  +  2xch4  +  Ufxfc )  ■  ^h2o(T)  +  xN2  ■  Iin2(T)  I  /( UfXfc ) 
'(1  -Ua)fio2(T)  +  3.7619fiN2(T)]  /(2 Ua) 

XH2hn2(.T)  +  XC0 JlcoCn  +  XcH4flCH4(T)  +  XCo2 hC02 (T)  +  XN2flN2(r)  +  (2xco  +  4xCH4)JlH2o(r)l  /(UfXfc) 
ho2(r)  +  3.7619hN2(T)]  /(2 Ua) 


however,  the  hybrid  efficiency  is  found  to  be  relatively  insensitive 
to  the  fuel  utilization  factor  and  the  SOFC  temperature  gradient. 
On  the  other  hand,  a  reduction  of  hybrid  efficiency  is  found  to  be 
evident  for  a  decrease  in  the  GT  temperature  ratio.  Other  factors 
affecting  the  hybrid  system  has  been  further  simulated,  such  as 
the  electrical  parasitic  consumption,  which  also  strongly  affects 
the  performance  of  the  combined  SOFC-GT  system.  For  instance, 
an  increase  of  9%  in  the  hybrid  efficiency  is  predicted  when  the 
power  consumption  factor  of  the  air  blower  decreases  from  20% 
to  0%,  indicating  that  effort  should  be  put  into  the  development  of 
more  efficient  blowers  to  reduce  the  system  parasitic  losses.  These 
interactions  analyses  can  be  used  to  determine  the  optimum  sys¬ 
tem  operation  and  thereby  to  derive  useful  design  guidelines.  It 
should  be  noted  that  the  size  of  system  is  not  defined  in  this  paper, 
because  the  model  developed  here  is  a  generic  one  which  is  able  to 
predict  the  performance  characteristics  of  a  wide  range  of  hybrid 
systems  potentially  sizing  from  2000  to  2500  Wm-2  with  efficien¬ 
cies  varying  roughly  between  50%  and  60%,  as  illustrated  in  the 
figures. 

In  summary,  the  hybrid  SOFC-GT  system  is  shown  to  be  robust 
and  thermodynamically  stable  over  a  wide  range  of  operating  con¬ 
ditions,  and  based  on  the  optimization  procedure  this  model  is 


2.  The  expressions  ofDi-Dio  in  Eq.  (56)  are  listed  as  follows: 

D  = _ MJe _ 

1  T)c  +  x(l  -  r)c)(T-[/T3) 

£>3  =  (1  -  *le)T3  +  D,Tr 

D5  =  1  +  Tq  -  D3 

D6  =  m,D4(D3-Ti  ^(g-l^)  +D2T3  (l  -  i)  (1  +  T0  -  D3) 
D7  =  im2  Ah(l  -  rife)  -  m3(T  -  T0) 
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Dg  —  1  —  TJe  + 


(l-r?c)D2(T1/r3)2 

T)e 


Dg 


(1  ~  '?c)D2(T1/T3) 

rie 


-Dj 


D10  =  T 


AT 

7T 


J3 

x 


Du 


1 

D^-Tq 


Dl2 


D5T3(1  -l/x)[(Tc-r3)/(D10x)-l] 

Tc-T3 


D i3  =  mj(Tj  - D3) 


(l-i)ln(D11)  +  D4D8D11 


Di4  =  miD4D8  ln(Dn)  +  D2  ^1  -  (D5  -  D8T3)  +  D12 

D15  =  -miD4(D9  +  l)ln(Dn)  -  miD4D9Dii(Ti  -D3) 

+  ^3  (1  -  d2d9 
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